Magnetic circular dichroism in GaMnAs: ( no ) evidence for an impurity band 
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Magneto-optical properties of the ferromagnetic semiconductor GaMnAs are studied in a material 
specific multi-band tight-binding approach. Two realistic models are compared: one has no impurity 
band while the other shows an impurity band for low Mn concentrations. The calculated magnetic 
circular dichroism (MCD) is positive for both models proving that, unlike previously asserted, the 
observed positive MCD signal is inconclusive as to the presence or absence of an impurity band in 
GaMnAs. The positive MCD is due to the antiferromagnetic p-d coupling and the transitions into 
the conduction band. 
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The ferromagnetic semiconductor GaMnAs [TJ |2j [3] 
is a fascinating material not only because of its impor- 
tance for technological applications [3], but also because 
its understanding poses great challenge. Indeed, under- 
standing the electronic structure of GaMnAs requires re- 
solving complexities due to disorder, correlations, and 
magnetism. Over the past years, an intense debate has 
sparked about the nature of the carrier states (holes) [5] : 
Do the relevant states reside in a GaAs-like valence band 
or in an Mn-induced impurity band? While the leading 
mean-field picture [2] of GaMnAs relies on the valence- 
band scenario, also supported by transport experiments 
[6] [7], recent optical experiments appear to imply the 
existence of an impurity band |5]. 

Ferromagnetism in GaMnAs is mediated by free holes 
induced by Mn impurities [9j [10]. Experimentally, 
magneto-optical studies of GaMnAs [8} HI} fT2 | [T3 l [T4l [T5] 
gave important insights into the exchange mechanism be- 
tween the holes and the local Mn moments [8]. In par- 
ticular, magnetic circular dichroism (MCD), which is the 
relative absorption strength for the left and right circu- 
larly polarized light, can give the signs and magnitudes 
of the exchange coupling constants p~6j [17]. It is now 
generally believed that the coupling between holes and 
the local moments (p-d coupling) is antiferromagnetic. 

The observed MCD signal around the fundamental ab- 
sorption edge, £? gap « 1.5 eV, is positive, which seems to 
imply a ferromagnetic p-d coupling, in contrast to what 
is known about II-Mn-VI diluted magnetic semiconduc- 
tors [TTJ [18] , as well as about GaMnAs at very low dop- 
ing levels, at which the material is paramagnetic [19]. 
Keeping the coupling antiferromagnetic, the natural ex- 
planation, within the valence-band picture, appeared to 
be a shift of the Fermi energy due to the Mn doping — Mn 
acts as an acceptor — in the ferromagnetic regime. This 
so called Moss-Burstein shift would then reverse the or- 
der in which the two circular polarizations of light are 
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On the other hand, recent MCD experiments have been 
interpreted to imply the existence of an impurity band 
[15], invalidating the Moss-Burstein shift picture. Such a 



view appears consistent with the calculations of Ref. [T8J 
which concluded that the positive MCD would arise in 
the valence-band scenario only if the p-d coupling were 
ferromagnetic (this argument counters the results of Ref. 
[22] which finds instead a positive MCD in the valence- 
band model with Mn disorder in the Born approxima- 
tion). Also, recent tight-binding calculations for an em- 
bedded single Mn impurity in a GaAs matrix [23 point 
to the existence of an impurity band. The apparent pic- 
ture coming from experiments [8] [131 H3 El] is that of a 
dominant and spectrally broad positive MCD signal from 
the transitions to an impurity band. 

Here we argue that the positive MCD signal has little 
to do with the presence or absence of an impurity band in 
GaMnAs. The signal results from the spin-resolved elec- 
tronic level ordering in GaMnAs which is present despite 
the strong disorder. The Fermi level shift due to dop- 
ing by Mn acceptors is strong enough to make the MCD 
signal positive, regardless whether or not the Fermi level 
lies in the valence band (Moss-Burstein shift) or in the 
impurity band, while holding to the antiferromagnetic p- 
d coupling in both cases. The dominant transitions are 
those involving the conduction band, not the impurity 
one. We thus find most previous conclusions drawn from 
the positive MCD signal unfounded. Our argument is 
based on large-scale tight-binding calculations of MCD in 
disordered GaMnAs using two models. One model repro- 
duces the valence-band picture, the other shows an im- 
purity band. Both have antiferromagnetic p-d coupling, 
and both show a positive MCD around the fundamental 
absorption edge, in agreement with experiment. 

Our simulations use a material-specific microscopic 
tight-binding approach. The models are based on 16 sp 3 
valence and conduction bands of GaAs which are approx- 
imated throughout the entire Brillouin zone to fit the ex- 
perimentally determined band structure [25] [26] . We use 
two different parameter sets for the inclusion of Mn im- 
purities into the models. One of the models [27] is char- 
acterized by an only slightly changed host valence band 
and a shift of the Fermi energy according to the number 
of holes that are added with the Mn. The other model 
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[28] shows a strong shift of the host valence-band states 
into the gap leading to the formation of an impurity band 
which starts to merge with the valence band at Mn con- 
centrations around 1% [29 . We present results on the ab- 
sorption of left and right circularly polarized light for fer- 
romagnetic bulk Gai-^Mn^As at zero temperature and 
various Mn concentrations. Within this framework we 
can treat the disorder effects non-perturbatively. While 
this approximation seems justified for large concentra- 
tions of Mn impurities [5 j, x > 1%, an explicit inclusion 
of the carrier-carrier interactions (computationally infea- 
sible) may give quantitative corrections of our findings 

As already stated, we investigate two different sets of 
tight-binding parameterizations to study in detail how an 
impurity band affects the MCD results. These two mod- 
els were thoroughly analyzed in Ref. [29] concerning the 
density of states, the size of the band gap, the position of 
the Fermi energy, localization properties of holes, optical 
effective masses, and mean free paths. The first model, 
which we call model A, was suggested by Masek [27] . 
This model, derived from a first-principles approach, de- 
scribes GaMnAs qualitatively very similar to what is ex- 
pected from a p-doped GaAs [29] . Its main characteristic 
is the inclusion of additional holes and therefore a corre- 
sponding shift of the Fermi energy into the host valence 
band. There is no formation of an impurity band within 
this model. 

The second model, which we call model B, was sug- 
gested by Tang and Flatte [28]. Within this model the 
Mn impurities are described by a modified on-site poten- 
tial and a spin-dependent potential at the four nearest 
As neighbors. The two relevant parameters are adjusted 
such that the experimental binding energy of 0.11 eV of 
the Mn impurity is recovered. Qualitatively, this model 
not only shifts the Fermi energy due to the additional 
holes but it also affects the positions of the energy lev- 
els of the host material rather strongly [29]. Applying 
this model to disordered systems with many Mn impuri- 
ties leads to the formation of an isolated impurity band 
for Mn concentrations x < 1%. At higher concentra- 
tions the impurity band starts to merge with the host 
valence band. Hence, the comparison of these two dif- 
ferent parameterizations, models A and B, allows us to 
draw conclusions on how the existence of an impurity 
band influences the MCD signals that are seen in the 
experiments. 

We analyze the MCD signal by evaluating the diagonal 
and off-diagonal elements of the optical absorption ma- 
trix. The matrix elements <j a p can be written in terms 
of the eigenstates \n) and eigenenergies E n as [31 
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FIG. 1: (Color online) Left panel: model A, right panel, 
model B. Upper panel: Schematic level ordering deduced from 
absorption data. Lower panel: Absorption of right (solid line) 
and left circularly polarized (dashed line) light. The system 
size is 2000 atoms with one Ga atom replaced by Mn. The 
absorption rates are evaluated at the T point in the super- 
lattice Brillouin zone. Levels are labeled according to their 
carrier character in pure GaAs: conduction band (CB), heavy 
and light-hole bands (HH, LH), and the spin-orbit split-off 
band (SO). Magnetic quantum numbers are indicated. Since 
substitutional Mn is an acceptor, the highest valence band 
level ('HH -3/2') is not occupied so that the transitions 'a' 
are suppressed. 



f n (Ep), energy difference A n / n = E n > — E ni and momen- 
tum matrix elements p^, n = (n f \p a \n). The absorption 
<j± for right (left) polarized light can be obtained by re- 
placing p — > p ± = p x ± \p y in Eq. ([!]). The MCD signal, 
for a constant index of refraction (or poor conductors), 
is 



MCD = (<T_ - Cr + )/((J_ + <7+). 



(2) 



with the system volume ^, the Fermi function f n 



In order to numerically evaluate the conductivity ([!]) 
we calculate the eigenenergies and the matrix elements 
of the momentum operator using a multi-band tight- 
binding approach [29 . This approach is applied to peri- 
odically repeated finite size super-cells containing up to 
2000 atoms. The resulting conductivity is furthermore 
averaged over several disorder configurations. Therefore, 
our approach goes beyond the single impurity calcula- 
tion presented on model B in Ref. [23] as we explicitly 
include disorder averaging effects. In our calculations, 
for numerical reasons, we use a smearing temperature 
corresponding to 1 meV in the Fermi functions and a 
peak width S = 5 meV. The integration over the super- 
lattice Brillouin zone is performed by summation over 
up to 2048 different k vectors [29]. As the systems un- 
der consideration are disordered one cannot restrict this 
summation to the irreducible part of the Brillouin zone. 
Furthermore, the finiteness of the systems limits the Mn 
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concentrations to x > 1% in our present simulations. 

For a better understanding of the numerical MCD re- 
sults we first analyze the effect of a single Mn impurity in 
a super-cell of 2000 atoms with periodic boundary condi- 
tions. The absorption was evaluated at the T point, i.e. 
k = in the super-lattice Brillouin zone. Due to the mag- 
netic moment of the Mn impurity the spin degeneracies of 
the conduction s- and the valence p-bands are lifted. The 
order in which the levels appear can be concluded from 
the absorption peaks of <j± at A « E gap ^1.5 eV. We 
show the numerical results in Fig. [I] The Fermi energy 
was chosen to be (artificially) above the highest impurity 
level, to show all relevant optical transitions. The ab- 
sorption peaks then correspond to the six possible tran- 
sitions between the six heavy hole, light hole, and spin 
orbit split-off bands to the two conduction bands [T7] . 
Other transitions are suppressed because of the selection 
rules for the momentum matrix elements. Among the al- 
lowed transitions, 'a' and 'c' are the strongest ones due to 
their larger momentum matrix elements. From the posi- 
tion of each <j + and <7_ peak one can uniquely identify 
the involved bands. The extracted ordering of the levels 
is shown schematically in the upper panel of Fig. [I] 

As a Mn impurity also acts as an acceptor, the high- 
est impurity level, that is the levels labeled 'HH -3/2' in 
Fig. [I] is in fact unoccupied. This means that the ac- 
tual Fermi energy lies just below this highest impurity 
level. Therefore, all transitions starting from this level 
are suppressed. For the absorption this means that the 
peak labeled 'a' is not observable at zero temperature. 

Figure [I] points to common features and differences be- 
tween the two models. The order of the first four peaks, 
'a' through 'd', is the same for both models implying the 
same ordering of the valence band states. The spin of 
the Mn core 3d electrons in our simulation points into 
the — z direction, that is, anti-parallel to the propagation 
direction of the light. This leads to an acceptor state 
with spin-down character implying antiferromagnetic p- 
d exchange coupling for the holes, for both models. 

Besides this common feature there are two major dif- 
ferences between model A and model B. First, the magni- 
tude of the level splitting due to the magnetic impurity is 
very different. For example, model A results in a splitting 
of the four heavy and light hole bands of « 4 meV while 
in model B this splitting is ~ 0.13 eV. This difference 
in the splittings is associated with the fact that model B 
explicitly reproduces the bound state energies of a single 
Mn impurity at « 0.11 eV in the host gap by a strong 
deformation of the host valence band [28] . The second 
major difference between the two models is the different 
order of the two conduction bands. In model A the spin- 
up band has a higher energy while in model B it is the 
spin-down band. Hence we find that model A is based on 
a ferromagnetic s-d coupling, consistent with the preva- 
lent view, while model B produces antiferromagnetic s-d 
coupling. Concerning the absorption data this is reflected 



^ 10 





















model 


A 


V 


"a 8 
"a 6 








O 5 




R 












--- 


a 








F 




1 










7 








Q 




\ i 












, 


















1 










h 




2 4 







































-2 






EleVI 








,y 





0.5 

Q 

U 

S o 
























































MCD 


























IMCD 












































^ 10 


















a 


















model 


B 




7 S 8 

"a 6 

CM 

2 4 

i o 








+ 




O 
















--- 


a 












/ 














i 


Q 























b 






















-2 






oartia 


1 












































MCD 


























0.4 

Q 

U 

S o 















































































































L5 A[eV] 2 



FIG. 2: (Color online) Calculated absorption of right (solid 
line) and left (dashed line) circularly polarized light for model 
A (upper graph) and model B (lower graph). The accompa- 
nying thin lines in model B are the partial conductivities ex- 
cluding the transitions to the conduction band. The system 
size is 432 atoms with 2 % Ga atoms replaced by Mn. The ab- 
sorption rates are evaluated using 2048 k vectors (256 for the 
partial conductivities) in the entire Brillouin zone; disorder 
averaging was over 5 configurations, and a running average 
over 25 meV was applied. From the absorption the MCD 
signal and the integrated MCD (IMCD) are determined. The 
insets show the density of states with the position of the Fermi 
energy (model A: E F = -0.23 eV; model B: E F = 0.13 eV). 
For model B, the energy Eo marks the top of the impurity 
band, E = 0.6 eV. 



in the fact that at energies A « E gap + Aso ~ 1-8 eV 
for model A the cr_ peak comes first while for model B it 
is the cr + peak. This specific in the ordering of the con- 
duction band does not alter the calculated MCD signal 
qualitatively. 

How are the above effects of a single Mn impurity re- 
flected in the absorption of the left and right circularly 
polarized light in strongly disordered samples with many 
impurities? The simulation results for a Mn concentra- 
tion x = 2% are shown in Fig. [2] For this Mn concen- 
tration model B still shows an impurity band, which is 
attached to the host valence band, while model A does 
not exhibit an impurity band at all as can be seen in 
the insets of Fig. [2] The Fermi energy in model A lies 
within the host valence band whereas for model B it lies 
within the impurity band. For both models we find that 
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ordering due to disorder. Our results are shown in Fig. [3] 
for x = 2% (E F = E ). The MCD is negative as the 
<j + transitions are now allowed, dominating the signal in 
accord with the single- impurity picture in Fig. [I] 

In summary, we find that the positive MCD signal in 
GaMnAs (a) is inconclusive as to the presence or ab- 
sence of an impurity band, (b) implies antiferromagnetic 
p-d coupling also for the valence-band scenario, and (c) 
originates from the transitions to the conduction band. 

We thank D. Weiss and C. Back for helpful discussions. 
This work was supported by the DFG SFB 689. 
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FIG. 3: (Color online) Top panel: model A, lower panel: 
model B. Integrated MCD signal for different Mn concen- 
trations x = 2% (solid line), x = 5% (short dashed), and 
x = 10% (long dashed). The system size is 432 atoms; the 
absorption rates were evaluated using 256 k vectors, and dis- 
order averaged over 20 configurations. For model B we also 
show the integrated MCD (IMCD) signal for Mn concentra- 
tion x — 2% with the Fermi energy shifted to the top of the 
impurity band, Ep = E (dashed dotted line). 



the cr_ absorption sets in at lower energies A implying 
a positive MCD signal. This is understandable as the 
strong a + transitions at low A « E gap are suppressed as 
the highest valence and impurity band states with spin 
down character are not occupied. Comparing with recent 
experiments we find very good qualitative agreement for 
the MCD and IMCD signals [12j[13]. A precise quantita- 
tive agreement cannot be expected from a tight-binding 
model [32] . 

To test the hypothesis put forward in [15 , that the 
positive MCD is due to the transitions to the impurity 
band, we have calculated the partial conductivities by 
excluding the conduction band states in Eq. [I] for model 
B; see Fig. [2] These conductivities (and their differ- 
ences) are much smaller than the complete ones, giving 
no support to the hypothesis. Instead, the positive MCD 
around the absorption edge is due to the transitions to 
the conduction band. 

As in Ref. 13, we show the IMCD signal for various 
Mn concentrations in Fig. [3] With increasing Mn concen- 
tration the magnitude of the IMCD signal is decreasing 
as the total absorption becomes stronger. However, the 
shape of the signal is rather robust for all investigated 
Mn concentrations. In particular, we do not observe a 
change in the sign of the MCD signal within the consid- 
ered range of concentrations x from 2 % to 10 %. 

What would the MCD signal be for a fully compen- 
sated GaMnAs? It is particularly interesting to see what 
the MCD were if the Fermi level were at the top of the 
impurity band in model B, as this reflects the ordering 
of the impurity band levels and the deformation of the 
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